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[i] Existing paleoclimate data suggest a complex evolution of hydroclimate within the 
Indo-Pacific Warm Pool (IPWP) during the Holocene epoch. Here we introduce a new leaf 
wax isotope record from Sulawesi, Indonesia and compare proxy water isotope data with 
ocean-atmosphere general circulation model (OAGCM) simulations to identify 
mechanisms influencing Holocene IPWP hydroclimate. Modeling simulations suggest 
that orbital forcing causes heterogenous changes in precipitation across the IPWP 
on a seasonal basis that may account for the differences in time-evolution of the proxy 
data at respective sites. Both the proxies and simulations suggest that precipitation 
variability during the September-November (SON) season is important for hydroclimate 
in Borneo. The preeminence of the SON season suggests that a seasonally lagged 
relationship between the Indian monsoon and Indian Ocean Walker circulation influences 
IPWP hydroclimatic variability during the Holocene. 
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1. Introduction 

[ 2 ] Variability in the hydrological cycle within the Indo- 
Paciftc Warm Pool (IPWP) is an important component of 
the global climate system: perturbations in deep tropical 
convection reverberate through oceanic and atmospheric 
teleconnections and affect both tropical and extratropical 
climates [e.g., Ropelewski and Halpert, 1987]. Many uncer- 
tainties remain regarding the dynamics of IPWP hydroclimate, 
including the response of Indo-Pacific Walker circulation to 
both past and future global temperature change [Vecchi et al., 
2006; Karnauskas et al., 2009; DiNezio et al., 2010, 2011]. 
Paleoclimate investigations that use both proxy data and 
modeling simulations provide an opportunity to improve our 
understanding of the hydrological cycle in the IPWP: by 
studying the response of IPWP hydroclimate to large, known 
climate forcings in the past - such as changes in the Earth’s 
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orbit - it is possible to both assess our ability to correctly 
simulate tropical climates and identify mechanisms that 
govern IPWP climate response to external forcing. 

[ 3 ] Paleoclimate archives that record water isotopes are 
useful tracers for IPWP hydroclimate, as water isotopes in the 
deep tropics track the amount of regional rainfall as well as 
synoptic-scale convective behavior [. LeGrande and Schmidt, 
2009; Lewis et al., 2010; Pausata et al., 201 1; Tierney et al., 
2011]. To date, there are two speleothem-based records of 
the oxygen isotopic composition of precipitation within the 
IPWP: one from Borneo, Malaysia [ Partin et al., 2007] and 
another from Flores, Indonesia [Griffiths et al., 2009]. 
The two records depict markedly different trends in hydro- 
climate since the last deglaciation, suggesting regional het- 
erogeneity within the IPWP that begs explanation. Partin 
et al. [2007] and Griffiths et al. [2009, 2010a] speculated 
that a mixture of forcing mechanisms - the influence of 
Earth’s precession on deep convection in the IPWP, the 
position of the Inter-Tropical Convergence Zone (ITCZ) and 
rising sea levels (for the earliest Holocene) - might explain 
the observed patterns. Any of these mechanisms is plausible, 
but none have been tested against general circulation model 
experiments. 

[ 4 ] Here we add a third IPWP isotopic record from off- 
shore Sulawesi, Indonesia, a site that sits precisely between 
Borneo and Flores (Figure 1). We use a different isotopic 
archive - the hydrogen isotopic composition of leaf waxes 

deposited in a marine sediment core - to reconstruct 
central IPWP hydroclimate for the last 15,000 years. Using 
isotopic information from three locations in the IPWP, 
we compare the proxy data to simulations conducted 
with coupled ocean-atmosphere general circulation models 
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et al., 2007], sedimentary bD lvaY from offshore Sulawesi (this 
study) and speleothem calcite and inclusion b 18 0 from the 
island of Flores [Griffiths et al., 2009, 2010b]. 


(OAGCMs), including the mid-Holocene (6K) simulations 
conducted as part of the Paleoclimate Intercomparison Project 
Phase II (PMIP2) and Holocene timeslice simulations con- 
ducted with the fully coupled, isotope-enabled Goddard 
Institute for Space Studies (GISS) ModelE-R. Using the 
proxy and model data in tandem, we discuss the mechanisms 
influencing IPWP hydroclimate and water isotopes during the 
Holocene epoch. 

2. Methods 
2.1. Proxy Data 

[5] We measured the hydrogen isotopic composition of 
leaf waxes (bD lvux ) in marine sediment core BJ8-03-70GGC 


(3.566°S, 119.383°E, 482 meters water depth), located 12 km 
off the southwest shore of Sulawesi in the Makassar Strait 
(Figure 1). Thirteen radiocarbon dates constrain the depth- 
age model for core 70GGC (Table 1) which is also described 
in Linsley et al. [2010] and Gibbons [20 12], The age model is 
based on a linear interpolation of the depth-age relationship 
between calibrated dates with the additional assumption that 
the top of the core is equivalent to 0 yr BP. Core 70GGC was 
sampled every ~100 years downcore, with somewhat higher 
resolution during the middle Holocene (near 4 ka) and the 
Younger Dryas chronozone (12.7-11.5 ka). SD wax analyses 
were conducted as previously described in Tierney et al. 
[2010], Briefly, freeze-dried sediments were extracted with 
a solvent mixture of d i c h I o ro m cth an can c t h an o 1 (9:l,v/v) 
using an Accelerated Solvent Extractor. Leaf waxes were 
purified from the resulting total lipid extract via NH 2 column 
chromatography, methylated with methanol of a known 
isotopic composition, and then further purified via silica gel 
chromatography. The C30 fatty acid was analyzed in triplicate 
for its hydrogen isotopic composition via gas chromatography- 
isotope ratio monitoring-mass spectrometry (GC-IR-MS), 
using a Thermo Delta XL mass spectrometer at Brown 
University. H 2 standard gas calibrated to VSMOW was injec- 
ted three times before and after the sample lipid peaks as an 
internal standard. In addition, an external fatty acid methyl 
ester (FAME) standard of known isotopic composition was 
run every nine injections to monitor drift. Isotopic values 
were corrected for the added methyl group. Results presented 
here are triplicate means, and average triplicate standard error 
was 0.7%o. 

[6] bT) wax is an emerging proxy for the isotopic composi- 
tion of precipitation and has been widely applied in tropical 
settings [e.g., Schefufi et al., 2005; Tierney et al., 2008], 
including at a nearby core site from the Makassar Strait 
[Tierney et al., 2010]. However, changes in wax source, 
including changes in vegetation or biomes through time, can 
complicate its fidelity as a recorder of the isotopes of pre- 
cipitation [Sachse et al., 2012, and references therein]. With 
regards to the geographic source of the leaf waxes, given the 
humid (precipitation rates of 2-3 m/yr) tropical forest envi- 
ronment of southern Sulawesi, the proximity of the core site 
to the shore (12 km, Figure 1), and the high accumulation 
rates along the Sulawesi margin attributable to terrestrial 


Table 1. Raw and Calibrated Planktonic Foraminifera 14 C Data for Core 70GGC a 


Accession 

Number 

Species 

Depth 

(cm) 

14 C Year 

1 (7 

Median Year 
BP 

Lower la 

Upper lcr 

OS-45435 

G. ruber & G. sacculifer 

63.5 

1310 

30 

853 

901 

816 

OS-54152 

G. ruber & mixed planktonics 

104 

2750 

50 

2472 

2540 

2352 

OS-54153 

G. ruber & mixed planktonics 

176 

4240 

50 

4341 

4413 

4269 

OS-60747 

G. ruber & mixed planktonics 

201.5 

4850 

40 

5155 

5245 

5079 

OS-60750 

G. ruber & mixed planktonics 

241.5 

5410 

35 

5789 

5847 

5737 

OS-54154 

G. ruber & mixed planktonics 

280 

6160 

60 

6594 

6666 

6513 

OS-60737 

G. ruber & mixed planktonics 

320.5 

7190 

45 

7653 

7695 

7598 

OS-70826 

mixed planktonics 

330.5 

8080 

40 

8538 

8589 

8481 

OS-70827 

mixed planktonics 

340.5 

8700 

45 

9380 

9438 

9325 

OS-70828 

mixed planktonics 

350.5 

9210 

45 

10030 

10124 

9954 

OS-60748 

G. ruber & mixed planktonics 

361.5 

10250 

40 

11223 

11251 

11187 

OS-65581 

mixed planktonics 

393.5 

12550 

50 

13990 

14064 

13898 

OS-45436 

G. ruber & G. sacculifer 

405.5 

12900 

45 

14649 

14934 

14257 


“From Linsley et al. [2010] and Gibbons [2012], The 14 C ages were calibrated with the Marine09 curve using Calib 6.0 [Reimer et al., 2009] with AR = 0. 
Depth-age pairs were linearly interpolated to every cm between dated intervals to form an age model, with the additional assumption that the coretop 
corresponds to 0 yr BP. 
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Figure 2. Observed and simulated seasonal cycles of precipitation for the Borneo, Sulawesi, and Flores 
sites (Figure 1). Blue lines represent the observed mean seasonal cycle of precipitation from the Global 
Precipitation Climatology Centre (GPCC) v5 gridded precipitation product (http://gpcc.dwd.de) [Rudolf 
et al., 2011], shading indicates the 68% confidence intervals. Black lines represent GISS ModelE-R 
simulated precipitation from the control experiment (OK). Red lines represents the simulated multimodel 
mean of 1 1 PMIP2 OAGCMs. For the simulated data, values are the mean of gridcells spanning 2-6°N, 
112.5-1 17.5°E (for northern Borneo), 2-6°S, 1 17.5-122.5°E (for Sulawesi) and 10°S, 1 17.5-122.5°E 
(for Flores). 


input (40-100 cm/ka [Oppo et al., 2009; Linsley et al., 
2010]) we assume that the waxes at Site 70GGC have 
always been locally derived from southern Sulawesi and 
predominantly transported by fluvial processes. Vegetation 
changes through time are a greater concern as pollen evi- 
dence from northern Sulawesi indicates expansion of grass- 
lands during the deglacial period [Dam et al, 2001]. Based 
on the apparent fractionation of studied grasses [Sachse et al. , 
2012], the effect of having grasses contribute proportionally 
more to the £>D wax signal during the deglaciation would be to 
attenuate the magnitude of the primary water isotope signal. 
From this perspective, the deglacial changes in 6D wax can be 
considered a minimum estimate. At around 10 ka, a forested 
woodland similar to today’s became established on Sulawesi 
[Dam et al., 2001] and lower-resolution analyses of the car- 
bon isotopic composition of leaf waxes in 70GGC indicate 
variability of less than 3%o (J. Saenz et al., manuscript in 
preparation, 2012) indicating a relatively stable Holocene 
landscape. These data suggest that SD wax can be confidently 
interpreted as a primary recorder of the <5D of precipitation 
across the timeframe of interest here. Nevertheless, our 
interpretation of the climatic meaning of the 5D wax data is 
deliberately earned out through tandem use of the existing 
speleothem and foraminiferal b 18 0 data and isotope-enabled 
modeling experiments. 

2.2. Model Simulations 

[ 7 ] We obtained simulated precipitation fields from the 
pre-industrial control (OK) and mid-Holocene (6K) experi- 
ments conducted as part of the Paleoclimate Intercomparison 
Project 2 (PMIP2) [Braconnot et al., 2007] from all available 
fully coupled OAGCM models archived in the PM1P2 data- 
base, which includes CSIRO-MkLl.l, 1PSL-CM4, CCSM3.0, 
FGOALS-gl.O, MIROC3.2, FOAM, MRI-CGCM2.3.4fa, 
MR1-CGCM2 . 3 .4nfa, GISS ModelE-R, ECHAM5-MPI, and 
HadCM3M2. The mid-Holocene experiments were forced 
with orbitally induced changes in solar insolation as well as 


greenhouse gases estimated from ice core data. Details of the 
models and the parameters used for the PM1P2 modeling 
experiments are available online at http://pmip2.lsce.ipsl.fr. 
We also use precipitation and hydrogen and oxygen isotopes 
of precipitation fields from time slice simulations conducted 
with one of the PMIP2 models, the isotope-enabled GISS 
ModelE-R. Time slice simulations were conducted for 9K, and 
then every IK from 6-OK. Similar to the PMIP2 protocol, 
these experiments were forced with insolation changes and 
greenhouse gases. The 9K GISS time slice simulation incor- 
porates a remnant Laurentide ice sheet, with mean ocean 
salinity and b 18 0 adjusted to account for the associated ice 
volume changes (35 psu, +0.33%o, equivalent to 35 meters of 
sea level change). Further details of the model, parameters 
used for the time slice simulations, as well as global results are 
described in LeGrande and Schmidt [2009]. 

[8] Though the sites of interest to this study (Borneo, 
Sulawesi, and Flores) are geographically close to one another 
(Figure 1), they experience different seasonal cycles in 
precipitation. GISS ModelE-R and the multimodel mean, 
respectively, reasonably approximate the distinct seasonal 
cycle of precipitation at these locales (Figure 2). This sug- 
gests that in spite of the relatively low spatial resolution of 
the PMIP2 OAGCMs, we can use the results of these 
simulations to investigate hydroclimate in this region with 
some confidence. 

3. Results and Discussion 

[ 9 ] The SD wax data from core 70GGC are relatively enri- 
ched (— 155%o) near 15,000 yr BP, then trend toward more 
depleted values (— 170%o) until ca. 8,000 yr BP, after which 
they stabilize (Figure 3a). Higher-frequency variability 
punctuates the 5D wax record throughout, including an 
enrichment event at ~900 yr BP during the Medieval Climate 
Anomaly (MCA) that replicates the result of Tierney et al. 
[2010], The data also show abrupt depletion around 11 ka, 
sustained depletion near 4 ka, and an enrichment event at 3 ka 
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0 5,000 10,000 15,000 

YEAR BP 

Figure 3. (a) Leaf wax hydrogen isotope data {b'D wax ) from 
marine sediment core 70GGC. Gray error bars represent the 
analytical 95% confidence intervals. Red line represents a 
500-yr gaussian smoothing of the data to emphasize long- 
term trends. Red triangles indicate the stratigraphic locations 
of 14 C dates (Table 1). MCA marks the Medieval Climate 
Anomaly and an enrichment event that replicates the result 
of Tierney et al. [2010]. Dotted lines mark the beginning 
and end of the Younger Dryas chronozone (YD), when no 
change in Sulawesi bD wax is observed. Other notable enrich- 
ment and depletion events as described in the main text are 
indicated, (b) Makassar Strait b 18 O sw , reconstructed from 
paired foraminiferal Mg/Ca and 6 1 8 0 measurements in several 
marine sediment cores including 70GGC [Lins ley et al., 
2010]. Gray error bars represent the compounded analytical 
and core-compositing uncertainty. 


(Figure 3a). There is no anomaly in SD wax associated with 
the Younger Dryas (YD) event, similar to the speleothem 
6 ls O data from Borneo [Partin et al., 2007], but unlike the 
speleothem <5 1 x O record from Flores, which show isotopic 
depletion during the YD [Griffiths et al., 2009]. 


[ 10 ] Comparison between bD wax and a reconstruction of 
seawater 6 18 0 (d 18 O i?w ) from the Makassar Strait (Figure 3b) 
indicates that both isotopic proxy records evolve similarly, 
with a linear trend toward isotopic depletion from 15-8 ka 
followed by relative stability during the rest of the Holocene. 
The 15-8 ka trend integrates several influences on 1PWP 
water isotopic composition, including changes in global ice 
volume, regional temperature, oceanic and atmospheric 
advection of fresh water into and out of the IPWP, and local 
rainfall amount. SST proxy data suggest that regional tem- 
peratures rose ca. 1.5°C across this time interval [Linsley 
et al., 2010], which would only induce a minor (ca. 1.5%o 
for bD, 0.19%o for 6 18 0) [Majoube, 1971] change in the 
isotopic composition of rainfall, and ice volume can only 
account for an additional ca. 4%o and 0.5%o for <5D and b 18 0, 
respectively (scaling inferred sea level changes to a glacial 
ocean enrichment of l%o) [Schrag et al., 1996]. Thus, the 
remaining ca. 9%o (0.7%o) change between 14 and 8 ka in 
the Sulawesi dD M , ax (6 I8 0. 5W ) record represents regional 
hydrological change. Seawater enrichment in the IPWP dur- 
ing the deglaciation and early Holocene is a robust feature of 
<5 18 0, w proxy records from the region [Stott et al., 2004; 
Oppo et al., 2007; Linsley et al., 2010] and has been 
explained by increased export of water vapor to the Asian 
continent and Indian Ocean as a result of an intensified 
monsoon, an increase in the advection of salty waters from 
the central Pacific, and reduced transport of water vapor across 
the Panama isthmus [Oppo et al., 2007; LeGrande and 
Schmidt, 2011]. The isotopic enrichment may also reflect a 
reduction in regional rainfall amount. Previous GISS ModelE- 
R isotope-enabled simulations [Lewis et al., 2010] identified 
rainfall amount as the climate parameter best correlated with 
precipitation isotopic variability in the Indonesian archipel- 
ago, and the timeslice simulations from GISS ModelE-R 
used here indicate that mean annual changes in precipitation 
amount and precipitation isotopes over Sulawesi during the 
Holocene (excluding 9K due to the simulated impact of ice 
volume) are strongly correlated (r = 0.94). Thus, the 15-8 ka 
trend in the Sulawesi SD wax data in particular likely reflects a 
shift from a drier deglacial period to wetter conditions during 
the Holocene. The simplest explanation for a dry deglacial 
period is a thermodynamic one (cooler temperatures cause 
reduced humidity and precipitation in the wet tropics) [e.g., 
Clement et al., 2004], but exposure of the Sunda Shelf may 
have promoted more severe regional divergence and drying 
[DiNezio et al., 2011]. 

[n] The absence of a trend from 8-0 ka in Sulawesi 8T) wax 
suggests that changes in the seasonal distribution of insola- 
tion (i.e., orbital precession) had little apparent effect on the 
isotopes of precipitation at this location and that meridional 
movement of the ITCZ - oft-invoked to explain rainfall 
isotopic evolution elsewhere in the tropics on this timescale 
[Burns, 2011] - is not relevant for this sector of the deep 
tropics. Likewise, speleothem b 18 0 records from Borneo and 
Flores depict a complex time-evolution of IPWP water iso- 
topes not immediately evocative of precessional influence 
(Figure 4). Like Sulawesi 6T) wax , these records show a trend 
from 15-8 ka, which suggests that deglacial forcings 
(changing ice volume, exposure of the Sunda Shelf) affected 
water isotopes uniformly across the Indonesian region. Yet 
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Figure 4. Simulated and proxy-reconstructed isotopic data 
from the IPWP ( Partin et al. [2007], Griffiths et al. [2009, 
2010b] and this study). Calcite <5 18 0 and 6T> wca data are plot- 
ted as 100-year gaussian-smoothed series to emphasize 
lower frequency variability, with analytical 95% confidence 
intervals. Fluid inclusion data from Flores [Griffiths et al., 
2010b] are plotted with lcr analytical error bars. Data are 
uncorrected for ice volume effects to facilitate comparison 
with the GISS ModelE-R timeslice results. Mean annual 
(solid lines) and SON (dotted lines) simulated precipitation 
isotopes from the GISS ModelE-R timeslice experiments 
are plotted on top of the proxy data; error bars indicate the 
lcr variation in the decadal means. For simulated isotopic 
values we show the precipitation-weighted mean <5 18 O or 
hD value from the same spatial boxes defined in Figure 2, 
normalized to the Holocene proxy mean for visual clarity. 
The anomalies in the proxy and model data are comparable. 
Dotted lines denote the Younger Dryas (YD) as in Figure 3. 


after 8 ka, the records diverge. At Borneo, the monotonic 
trend in b 18 0 continues until 4 ka with a change in h 18 0 
between 10-4 ka of ~1.5%o (Figure 4). In contrast, the h 18 0 
record from Flores shows millennial-scale oscillations during 
the Holocene and more enriched <5 18 0 near 4 ka (Figure 4). 

[ 12 ] In addition to calcite <5 ls O, Griffiths et al. [2010b] 
measured the b 18 0 of fluid inclusions within the Flores spe- 
leothem, which is presumably representative of the isotopic 
composition of the source water from which the calcite pre- 
cipitated (Figure 4). While fluid inclusion h l8 0 shares the 
deglacial trend with the calcite d 18 0, its evolution is different 
during the Holocene in that it indicates more depleted values 
at 4-5 ka. The difference between calcite b 18 0 and inclusion 
b 18 0 may reflect the effect of changing cave temperatures on 
the fractionation between drip water and calcite. During the 
mid-Holocene, this difference is 0.7%o, which would theo- 
retically imply that temperatures in the Flores cave were 
~3°C colder than present [Griffiths et al., 2010b]. This pos- 
sible temperature effect is significant in the context of the l%o 
variability in Holocene calcite b 18 0. (Figure 4). Although 
the errors on the fluid inclusion measurements are large 
(Figure 4) and thus allow for a wide range in possible cave 
temperature estimates, these data present a possibility that 
cave temperatures influenced Flores calcite <5 18 0 signifi- 
cantly, independent of rainfall b 18 0. Given this possibility 
along with the lack of replication of the record in a different 
cave, our conservative interpretation of the b 18 0 data from the 
Flores is that they do not resolve a mean change in local 
rainfall <3 1 8 0 across the Holocene. In this sense, the Flores 
data are in rough agreement with the SD wax data from Sula- 
wesi. While cave temperatures could have influenced the 
Borneo record as well, unlike Flores the Borneo data were 
replicated in different caves 5 km apart [Partin et al., 2007]. 
A temperature overprint would thus require coordinated 
temperature variability between caves, presumably dictated 
by large changes in air temperatures. As SST proxy data 
suggest that temperatures were relatively stable in the IPWP 
across the Holocene [Lins ley etal, 2010], this seems unlikely. 
In sum, the Holocene IPWP proxy isotopic data suggest peak 
depletion in b 18 0 at Borneo from 4-5 ka, but no detectable 
change in the mean isotopic composition of rainfall at 
Sulawesi and Flores. 

[ 13 ] If we accept that the differences between the proxy 
data sites are correct, the question then arises, why are the 
water isotope data different between these sites? We turn to 
simulations of IPWP hydroclimate in coupled OAGCMs to 
seek answers. First, we use the GISS ModelE-R timeslice 
simulations to directly compare simulated water isotopes with 
the proxy data (Figure 4). Weighted mean annual precipita- 
tion isotopes are enriched at 9K relative to OK but then change 
little across the rest of the Holocene at all three sites, in gen- 
eral agreement with the proxy data at Sulawesi and Flores, but 
not at Borneo. Rather, calcite h 18 0 at Borneo shows a better 
match with the simulated SON-season oxygen isotopic com- 
position of precipitation (h 18 O p ), both of which show maxi- 
mum depletion at 4 ka (Figure 4). An observational study of 
the seasonal cycle of isotopes in precipitation in northern 
Borneo found that the time of maximum depletion in precip- 
itation isotopes occurs during SON, when much of the water 
vapor originates from more remote southeasterly sources 
[Cobb et al., 2007]. Thus, rainfall variability during SON can 
strongly influence weighted mean-annual h^O^. GISS 
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Figure 5. Comparison between an ad-hoc estimation of 
the seasonal cycle in precipitation b 18 0 (6 1 8 Op) in northern 
Borneo, based on a three-year time series of observed rainwater 
and dripwater b 18 0 (black) [from Cobb et al., 2007] and the 
simulated seasonal cycle in <5 18 0p in northern Borneo by 
GISS ModelE-R (blue). 

ModelE-R simulates this seasonal cycle in () 18 Op at Borneo 
remarkably well (Figure 5), and as is the case with the 
observational data [Cobb et al., 2007], simulated <5 18 0p has 


only a weak correlation ( r = 0.28) with the seasonal cycle in 
rainfall amount (in contrast, simulated precipitation and 
b 18 O p are strongly correlated at both the Sulawesi and Flores 
sites, r = 0.84 and r = 0.92, respectively). The mismatch 
between simulated mean annual 6 ls O p and proxy <5 ls O at 
Borneo, then, does not reflect a fundamental deficiency in the 
model, but rather may simply indicate that the model is 
underestimating the amount of SON rainfall at 4K and/or the 
contribution of water vapor from the southeasterly source. 
Regardless, the depletion at 4 ka (Figure 4) seen only in the 
Borneo proxy record suggests a hydroclimatic response that 
is seasonal in nature and primarily affects the Borneo land- 
mass and immediate vicinity. 

[ 14 ] What mechanisms might cause a unique hydroclimatic 
response in Borneo? Though we do not have multimodel 
simulations specifically at 4 ka, the multimodel mean of 
6K precipitation simulated by 1 1 OAGCMs participating in 
PMIP2 (which includes GISS ModelE-R) yields insight into 
robust seasonal changes in mid-Holocene 1PWP climate that 
may partially explain the heterogeneity in the proxy data. On 
a mean annual basis, the simulated changes in precipitation at 
6K relative to OK are close to zero (+0. 1 mm/day, 6K-0K), in 
general agreement with the proxy isotopic data specifically at 
6K (Figure 4). However, on a seasonal basis there is a 
heterogenous response of hydroclimate to mid-Holocene 
climate conditions between Borneo, Sulawesi, and Flores 
(Figure 6). In boreal summer (June-July- August, JJA), pre- 
cipitation increases at 6K relative to OK over large landmasses 
in central Indonesia, causing wetter conditions over Borneo 
and Sumatra. In contrast, Sulawesi and Flores experience little 
change. Conversely, in boreal winter (December-January- 
February, DJF) precipitation decreases over all large land 


6K-0K JJA AP (mm/day) 6K-0K SON AP (mm/day) 





-1 - 0.8 - 0.6 - 0.4 - 0.2 0 0.2 0.4 0.6 0.8 1 


Figure 6. Multimodel mean of seasonal changes precipitation for the mid-Holocene (6K) relative to the 
preindustrial control (OK), as simulated by 11 PMIP2 OAGCMs: CSIRO-MkLl.l, IPSL-CM4, CCSM3.0, 
FGOALS-gl.O, MIROC3.2, FOAM, MRI-CGCM2.3.4fa, MRI-CGCM2.3.4nfa, GISS ModelE-R, 
ECHAM5-MPI, and HadCM3M2. B, S, and F = locations of the Borneo, Sulawesi, and Flores proxy 
isotope data. Stippling indicates regions where all 1 1 models agree on the sign of change. 
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Figure 7. The 9K versus 4K anomalies in precipitation and 
the oxygen isotopes of precipitation simulated by GISS 
ModelE-R. Stippling indicates significant (p < 0.05) changes 
based on comparison with average decadal variability. B, S, 
and F = locations of the Borneo, Sulawesi, and Flores proxy 
isotope data. 


surfaces in the IPWP - both north and south of the equator - 
and increases over the Indian ocean (Figure 6). Flores, as a 
small island, experiences slightly wetter conditions during 
this season as precipitation shifts from land to sea, whereas 
Borneo is slightly drier. Generally speaking, the seasonal 
redistribution of precipitation at 6K in the Indo-Pacific shows 
little meridional structure; rather, the changes in precipitation 
are from land to sea and vice versa. This reflects a differential 
influence of enhanced insolation seasonality, wherein 
increased insolation during the JJA season at 6K heats land- 
masses faster relative to the ocean, enhancing monsoonal 
circulation throughout the northern hemisphere and in the 
equatorial tropics, and decreased insolation during the DJF 
season at 6K cools land faster, resulting in divergence and 
enhanced export of moisture to areas over the ocean near the 
major landmasses [Liu etal, 2004; Clement etal., 2004; Zhao 
et al., 2005], 

[is] These model results suggest that Borneo’s unique 
hydroclimatic history could partially reflect the fact that the 
island is a large enough landmass to respond to temperature- 
driven changes in seasonal convection in-phase with the 
other landmasses in the IPWP. However, if the JJA season 
dominated annual-mean 6 n O p at Borneo (as it arguably did 


on mainland Asia) [ Wang et al., 2005], we might expect that 
Borneo b 1 8 0 would show peak depletion near 11-9 ka, 
when orbital configuration results in perihelion during JJA. 
Yet both the proxy <5 ls O and simulated <5 ls O from GISS 
ModelE-R suggest maximum depletion at 4K and not 9K. 
The 4 ka minimum may reflect a direct influence of 
increased insolation during the SON season (orbital config- 
uration puts perihelion near SON at 4 ka), as Partin et al. 
[2007] hypothesized. Here we explicitly test that hypothesis 
by comparing the 9K and 4K results from the GISS ModelE-R 
timeslice experiments. We find that indeed, precipitation 
amount is somewhat reduced over Borneo at 9K relative to 4K 
in the GISS ModelE-R simulations and that the isotopic 
composition of rainfall is enriched by ca. l%o (Figure 7). The 
latter anomaly is in excellent agreement with the proxy data 
(Figure 4), but is larger then expected from the amount effect 
alone, implying less rain derived from the remote southeast- 
erly source and possibly weaker convection at 9K relative to 
4K. However, there is no indication that the drying over 
Borneo at 9K relative to 4K during SON was caused by 
cooling associated with the change in seasonal insolation; 
there is a negligible anomaly in surface air temperature 
(ASAT 9K-4K = -0.02°C). 

[ 16 ] Alternatively, we argue that the timing of peak 
depletion in b 18 0 at Borneo at 4K can be explained in part 
by large-scale changes in Walker circulation across the 
Indian Ocean basin. In the 6K PMIP2 simulations, there is a 
zonal redistribution of precipitation across the Indian Ocean 
basin during SON, with an increase in precipitation over the 
western Indian Ocean and a decrease over the eastern Indian 
Ocean off of Sumatra (Figure 6). This anomaly and the 
mechanisms that cause it have been discussed in the litera- 
ture previously [Zhao et al., 2005; Tierney et al., 2011]. 
Tierney et al. [201 1] argued that the change in convection is 
a seasonally lagged effect of Indian monsoon circulation on 
the western Indian Ocean heat budget. In JJA, the Indian 
monsoon winds are shifted northward during the early to 
mid-Holocene, reducing latent heat export out of the western 
Indian Ocean. This feeds the development of a warm SST 
anomaly several months later, which triggers an atmospheric 
response resulting in increased convection and rainfall in the 
western Indian Ocean, at the expense of divergence and 
drying in the eastern Indian Ocean in the ascending branch 
of the Walker cell. The effect of this mechanism on the 
hydrological cycle in the Indian Ocean is large and robustly 
simulated by different models (Figure 6). The zonal pattern 
of anomalies in precipitation and the isotopes of precipita- 
tion at 9K relative to 4K in GISS ModelE-R suggests that 
this mechanism is stronger in the early Holocene than at 4 ka 
(Figure 7), as might be expected as it is related to the posi- 
tion and strength of the Indian monsoonal winds during 
boreal summer (which are stronger/farther north at 9K). 
Upper-level (250 mb) zonal westerly anomalies between 9K 
and 4K (0.9 m/s) and easterly zonal surface wind anomalies 
(—0.6 m/s) further suggest that there is a relative weakening 
in Indian Ocean Walker circulation at 9K during SON. 

[ 17 ] Given these results, we hypothesize that the observed 
trend in Borneo hydroclimate across the Holocene reflects a 
balance between two competing mechanisms: the direct 
effect that JJA perihelion has on precipitation due to heating 
and convection, and the indirect effect that JJA perihelion 
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has on coupled ocean-atmosphere dynamics during SON. In 
the early Holocene, the latter mechanism dominates, 
whereas in the mid-Holocene, the former mechanism dom- 
inates, resulting in maximum isotopic depletion at a point of 
compromise around 4 ka. Both of these mechanisms are less 
important for hydroclimate at Flores and Sulawesi, which do 
not experience increased precipitation during JJA (Figure 6), 
and - if the GISS results are a robust indication - sit outside 
the core of the isotopic anomaly associated with weaker 
Indian Walker circulation (Figure 7). 

4, Conclusions 

[is] Our investigation of model and proxy evidence for 
changes in 1PWP hydroclimate during the Holocene reveals 
that IPWP response to orbitally induced changes in insola- 
tion is heterogenous across the region on the seasonal level, 
and thus differences in hydroclimatic evolution between 
Borneo and central Indonesia are expected. The Sulawesi 
site in particular sits near the zero line of most of the sea- 
sonal anomalies (Figure 6) and so the new proxy 6D wax data 
we present here - which imply little change in mean bD of 
precipitation across the Holocene - are in good agreement 
with the model simulations. As temperature changes com- 
plicate the interpretation of Flores calcite b 18 0, we cannot 
robustly determine whether the model simulations agree 
with the proxy data. At Borneo, we find excellent agreement 
between the £ 18 0 proxy data and simulated SON 6 1 8 0 of 
precipitation in GISS ModelE-R, and the Borneo site is 
uniquely sensitive to large-scale changes in IPWP hydrology 
that occur during this season. The development of longer 
proxy isotopic records from Borneo and the vicinity will be 
an important test of our hypothesis that two competing 
mechanisms - both related to Indian Ocean atmospheric 
circulation - influence Borneo hydroclimate on precessional 
timescales. 
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